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We discuss ferroelectric phase transitions into single- and multidomain states in very thin films 
using continuous theory. It is shown that in nearly cubic ferroelectrics the domain state may survive 
down to atomic film thicknesses, unlike the single domain state, which is almost always unstable 
or metastable. This conclusion is valid almost irrespective of the nature of electrodes (metallic or 
semiconducting) and whether or not the screening carriers may be present in the ferroelectric itself. 



With a thrust towards developing nanoelectronics com- 
ponents, like ferroelectric (FE) memories [1] the question 
of phase transformation and the very existence of ferro- 
electricity (problem of "critical thickness" , first raised 
decades ago [2-5]) becomes ever more important. It is 
a focus of experimental and theoretical research, see e.g. 
[6,7]. Here we will present the thermodynamic results 
for phase behavior of ferroelectric thin films with various 
electrodes and consider both homogeneous and inhomo- 
geneous states of the film as a function of temperature 
and thickness. Many of prior works have given some re- 
sults almost exclusively for homogeneous (monodomain) 
ferroelectric film states [2-5,8] with estimates for a "crit- 
ical thickness" ranging from ~ 1/im [2] to a few nm [8]. 
We shall call critical the thickness where either homo- 
geneous or inhomogeneous spontaneous polarization be- 
comes unfavorable. Obviously, it is non-universal, and 
depends on many details like the ferroelectric material, 
lattice mismatch, type of electrodes, etc. For instance, we 
show below that the homogeneous state studied in Ref. [7] 
may exist only when very special conditions are met, in 
most real situations the FE films are likely metastable 
or even unstable with regards to breaking into domains. 
Importantly, ferroelectric domains in nearly unstrained 
cubic FE films may exist down to practically one unit 
cell. 

Consider a slab of uniaxial ferroelectric occupying the 
region —1/2 < z < 1/2 between two metallic (or degen- 
erate semiconductor) electrodes. The polar axis is per- 
pendicular to the film plane (parallel to the z— axis), and 
the bias voltage is set as —(+)U/2 at the right (left) 
electrode. The electric field would penetrate into elec- 
trode over very short Thomas- Fermi screening length [9] , 
and this screening layer manifest itself as some addi- 
tional "interface capacitance" [10]. To characterize the 
screening, we introduce the band bending potential 4> as 
ip = <p =F U/2. Since the screening is assumed to be very 
efficient, the potential drop in the electrode is small com- 
pared to the Fermi energy /i, \q<fi\/ /j, <C f , and the Poisson 
equation linearizes to <j>" = K 2 (j), where k 
is the inverse screening length, k~ 



y/6imq 2 /e s n 
0.5 — 0.7A in metals 



[9], comparable to the atomic distance d a t, n the net car- 
rier density, e s the dielectric constant, q the elementary 
charge. According to this equation, the band bending 
is antisymmetric, 4>{—z) = —cf>(z), </> = a^e""^ - '/ 2 ) at 
z > 1/2, with 



01 = (2ttPI + U/2) /(l + e s Kl/2). 



(1) 



This is the local electrodynamic boundary condition. In 
fact, the surface has properties different from the bulk 
and generally produces an effective "field" w coupled to 
the normal component of polarization, as was pointed 
out recently [11]. This, in principle, requires to add to 
(1) the additional boundary conditions for polarization 
[12,13,10,11]. The solution (1) applies when \qc/>i\ <C jit, 
that translates for metal, where nl ~ l/d a t 3> 1, into 
Anq\P\/n ~ 4irq\P\d at <C ft, which is always satisfied. 
The condition would be violated in moderately/lightly 
doped semiconductors. 

The equation of state of a ferroelectric 



AP + BP 3 = E, 



(2) 



where the homogeneous field in the monodomain ferro- 
electric would be E — (e s K,U - 8ttP) / {e s nl + 2) , A = 
(T—T c )/Tq, T c is the critical temperature, To the charac- 
teristic temperature, To ~ T at in displacive ferroelectrics, 
where T at ~ 10 4 — 10 5 K is the characteristic "atomic" 
temperature. The field E = E + Ed, where E = U/l is 
the external and 



E d = -PLo/l 



(3) 



is the depolarizing field, Lq — 87r/e s n the characteris- 
tic length scale in electrodes, and we have used the fact 
that in a metal the screening length is small, nl 3> 1. 
The homogeneous ferroelectric state in short-circuited 
electrodes, U = 0, requires an intersection of the P (E) 
curves given by Eqs. (2) and (3) at P ^ (Fig. 1). This 
necessary, but not sufficient, condition is realized at 



A < A h = -L /l, 



(4) 



illustrated in Fig. la, where curves 1-3 are plotted for 
Ti > T2 > T3. The nontrivial solutions exist only for 



1 




Ti>T 2 >T 3 
I = const 



E 



(b) 




\ 



const 



metastable 
unstable 



E 



FIG. 1. Solutions of the equation of state for FE film (2) 
with the depletion field Ed (3): (a) at constant thickness of 
the film I and decreasing temperature; (b) at constant temper- 
ature T and varying thickness I. Nontrivial solutions, when 
they exist, are either unstable (2) or metastable (3) with re- 
spect to domains. 



T2 : 3. If these solutions were stable with respect to forma- 
tion of inhomogeneous states, the corresponding critical 
temperature would be suppressed by AT C = —TqLq/1 ~ 
87rTod at /e s Z, compared to the bulk T c [4,10]. Hence, the 
homogeneous ferroelectricity would become impossible in 
films with a thickness below 



Ihc = 87r/e s Av|v4| max = LqTq/T c , 



(5) 



d at we find 



where |^4| max = T c /T . In displacive systems |A| 
1(T 2 - 1(T 3 , so for e s ~ 1 - 10 and n m ~ d~ t x 
l hc ~ (10 2 - 10 3 ) d at ~ 50 - 500A. One should bear in 
mind that the ferroelectricity was observed in thinner 
films [6] but it was not possible to check if the film was 
mono- or polydomain, the latter being much more likely 
scenario. Obviously, when the external field is larger than 
the depolarizing field (3), Eo > Ed, the film will be in 
monodomain state. This indicates the boundary between 
poly- and monodomain states and correctly describes the 
tilt of the hysteresis loops with decreasing film thickness 
observed in Ref. [6]. 

The necessity to study the sufficient conditions is sug- 
gested by Fig. 1, since the homogeneous ferroelectric so- 
lution belongs in the regions which correspond to either 



unstable or metastable states in the bulk samples. The 
problem of stability loss with respect to small fluctua- 
tions is a nontrivial yet tractable (linear) problem [14,15]. 
The check on metastability requires the calculation of free 
energies, which are much more difficult to find, so only 
limiting cases will be discussed below. The question of 
stability with respect to small fluctuations splits into two 
parts: stability with respect to (i) homogeneous and (ii) 
inhomogeneous fluctuations. Any solution of Eqs. (2) and 
(3) is stable with regards to homogeneous fluctuations. 
One can prove this by solving a relaxation dynamics for 
polarization constructed by generalizing the equation of 
state (2): 



BP 

7 - = -AP-BP 3 , 



(6) 



where the relaxation parameter 7 > 0. Linearizing this 
equation about the point P — Pq, where 



P =J-A/B 1 A = A + L /l 1 



(7) 



one obtains for homogeneous fluctuations 6 h P = P- P a 



7" 



dS h P 



ot 



A + 3BP 2 6 h P, 



(8) 



so that the perturbation does indeed decay with time as 
5 h P oc e at with the decrement a = 2A < 0. 

Now let us discuss stability with respect to inhomo- 
geneous fluctuations of polarization. The present situ- 
ation with the field penetrating into electrode over the 
Thomas-Fermi length is analogous to a ferroelectric film 
separated from the electrodes by nonferroelectric "dead" 
layers [16,17,14] that promote breaking the film into do- 
mains. To obtain the conditions for domain instability in 
the present case, one looks for a nontrivial solution of the 
equation of state with the gradient term included [15] 



AP + BP 3 - gV\P = E. 



(9) 



where V 2 i = <9 2 



d 2 
v 



is "in-plane" Laplacian, together 



with Maxwell equations. The solution is sought in a form 
of a "polarization wave" P = Pq + SPk(z)e lkx , <p — 
(pk(z)e lkx , where 5Pk(z) is the small perturbation over 
the homogeneous polarization Pq [14, 15]. In the case of 
metallic screening we obtain the following condition for 
an existence of the nontrivial solution with certain k : 



Xtan>= ^ 



(10) 



where x 2 = -e±(A + gk 2 )/4-K > 0, A = A + 3BP 2 . The 
case of interest to us is k > fc, which is easy to meet in 
the present case of metallic electrodes. We assume (and 
check validity later) that e±k/e s KX ^ 1- Then, the equa- 
tion simplifies to \kl = tt, the same as in FE film with- 
out electrodes. Substituting there the above expression 



2 



for x(k, A), we find the maximal value of A c = —2gk 2 at 
k = k c where this equality is first met, defining the tem- 
perature where the instability sets in. The "polarization 
wave" forms at A < such that 



linear cubic ~ 9 / ^ ~ (d a t / 1) 2 ■ 



(14) 



-A c 



2gk 2 c 



V±gl 2 ) 



A d = X/l, 

1/4 



(11) 

(12) 



where g is the coefficient before the gradient term in the 
equation of state (9), A = 4tt 3 / 2 5 1 / 2 /e[ 2 ~ d at le[ 2 
is the characteristic thickness for the domains forma- 
tion. Now, checking the assumption that we used 
to solve the Eq. (10), we see that it boils down to 

\/47reX/ (esy/Kpg'j ~ \/47te±/e s > 1, which is met in 

(e.g. perovskite) films with e± ~ 100 — 1000, and the 
typical e s ~ 1 — 10. 

These results remain basically unchanged if one were 
to account for realistic boundary conditions and the fact 
that the interface creates the effective "field" w coupled 
to the normal polarization component [11]. The surface 
field results in a "frozen" polarization Po(z), but since 
it extends over a short length on the order of a lattice 
spacing, it makes the film response "hard" just near the 
boundary. The film remains "soft" in the bulk and there 
the stability loss proceeds by the above scenario via ap- 
pearance of the polarization wave. The effect of the sur- 
face field is basically that the surface layer of atomic 
thickness would be excluded from the process of domain 
formation. 

Importantly, in nearly unstrained cubic ferroelectric 
films the multidomain states are much more likely than 
in standard uniaxial ferroelectrics. Indeed, upon cooling 
from the paraelectric phase of a uniaxial ferroelectric film 
a domain instability sets in a form of a "polarization 
wave" when A = Ad (11) with the critical wave vector 
k = k c (12) found from 



|A d | uni a Xial = 4^/y/ 2 / e f 



(13) 



Eq. (11), where ej_ — l + Air/A± is the dielectric constant 
in the direction perpendicular to ferroelectric axis in the 
plane of the film. This is the case of interest to us, since 
e_L ~ 10 3 in BaTi0 3 at 2% lattice misfit [18,7]. In stan- 
dard situations Eq.(13) applies, but if cubic perovskite 
films are grown on a substrate with small to negligible 
lattice mismatch, the period of a domain structure in a 
cubic ferroelectrics would increase. At the limit of appli- 
cability of this formula one can estimate A± w An <~ gk 2 
and we will call them "near-cubic" . Substituting this 
into (13), we obtain k c ~ 7r/(2 1//2 Z) ~ l/l, so that the do- 
main width a <~ Trfc^ 1 tends to become comparable to the 
thickness of the film. In this case, according to (13), the 
transition occurs very near the bulk critical temperature 
at 



The minimal thickness, where the near cubic thin 
films can still transform into a polydomain ferroelec- 
tric state, is then estimated for e.g. BaTiC-3 as Id <~ 
d at (max |A d | noar cubic)~ 1/2 ~ 6d at ~ 10A at low tem- 
perature. Therefore, the ferroelectricity in near cubic 
ferroelectric films can exist in a polydomain form down 
to "atomic" thicknesses (just about one unit cell thick), 
where the present continuous theory is at the border of 
validity but still able to produce semiquantitative results. 
One can estimate from (13) that this regime corresponds 
to e± w Al 2 /g ~ (l/d at ) 2 . A more careful consideration 
suggests that the transition in a cubic ferroelectric film 
with a smaller lattice misfit with the substrate than the 
above borderline value proceeds into a monodomain state 
with homogeneous in-plane polarization, e.g. P x . This 
threshold misfit is very small indeed, so in most cases the 
film remains uniaxial and splits into domains according 
to the above scenario. 

One can identify two possible cases from Eqs. (4) and 
(11), while considering possible phase transitions upon 
lowering temperature: (i) Ah > Ad, when A > Lq, the 
phase transition is into a homogeneous state, and (ii) 
Ah < Ad when A < Lq, the phase transition proceeds 
into a multidomain state. The case (i) does not mean 
that the state remains homogeneous. It is stable with re- 
spect to inhomogeneous fluctuations, but when the tem- 
perature approaches A w Ad, domains would start to 
form. Indeed, at this point the domain wall thickness 
W <~ \j 1 —gjAd becomes smaller to the domain width a, 
W < a, and the existence of usual domains becomes pos- 
sible. The situation with metallic electrode is analogous 
to a system with dead layers, which is metastable far from 
the critical point at \A\ <C \Ad\ with respect to domain 
structure at any thickness of the dead layer [16], and the 
domains are likely to form already at A w A d - Consider 
now the case (ii) A < L , where paraphase becomes un- 
stable with respect to domains close to the phase transi- 
tion point. Interestingly, it was found in Ref. [14] that the 
homogeneous state may be stable with respect to small 
inhomogeneous fluctuations not only in the paraphase 
but also in the ferrophase at some T <T S . This result is 
formal, however, since the system is likely metastable in 
this region and domains will grow. The temperature T s 
is defined by the condition A + SBPq = A d = —A//. 
For the case g = A = this corresponds to a point 
where a depolarizing field equals thermodynamic coer- 
cive field, separating bulk homogeneous metastable from 
unstable states. At T < T s and negligible energy of the 
domain walls (g = 0) the depolarizing field would split 
the FE film into domains. Therefore, the homogeneous 
state below T s is actually metastable, and this is also easy 
to prove by comparing its free energy with the domain 
state. The case with g ^ is more involved, but should 
be qualitatively similar, so that the low-temperature ho- 



3 



mogeneous state predicted in Ref. [14] is actually unob- 
servable in the thermodynamic sense. 

It has been found in first-principles modeling of 
BaTi03 ferroelectric film with SrRuC>3 metallic elec- 
trodes with in-plane lattice parameter corresponding to 
a SrTiC>3 substrate (i.e. mimicking a capacitor struc- 
ture grown on top of SrTiOa with ~ 2% compres- 
sive strain) that the state with a homogeneous polar- 
ization in c— direction remains stable down to ~ 24A 
at zero temperature [7] (~10A in unstrained PbTiC>3 
[19]). Further, it was claimed that the depolarizing field 
was solely responsible for vanishing of the ferroelectricity 
in thinner films, the chemistry of the interface (within 
the present nomenclature: the additional boundary con- 
ditions, ABC) was deemed unimportant (however, see 
[19]). 

With regards to such phase transitions with thickness 
of the film, there are again two possibilities for samples 
with the thickness I = lh = Lq/\A\ where Pq = 0: (i) 
paraelectric state with zero polarization is stable if A — 
Ah > Ad (i.e. L < A), and (ii) it is unstable and domains 
form, L > A. The condition (i) is met when 



< 1. 



(15) 



For very moderate values for perovskites e± ~ 100—1000, 
the electrode dielectric constant should then be e s > 
10 — 30, which may or may not hold for particular elec- 
trodes. This condition has not been checked in [7] and, 
therefore, it remains unclear what regime corresponds to 
the phase transitions studied in [7]. We see that the 
assumption about the existence of homogeneous state 
made in [7] is highly questionable and most likely the 
film in the ground state would be split into domains. 
When the condition (15) is met [case (i)], the film is in 
a paraelectric phase when I < lh (see, however, a reser- 
vation below) and becomes single domain ferroelectric in 
thicker films, I > lh, Fig. lb. Indeed, at I > lh we have 
A-A d = 2L a (jj- - +^f^ > and the homogeneous 
state remains stable with regards to small inhomogeneous 
fluctuations. However, it should become metastable at 
I > Id, where domains become possible, i.e. W < a, 

which happens at I > l d , Id ~ \[izi^ ~ > lfl - 
The paraelectric phase is stable at I < lh with respect 
to small fluctuations but may, at least in some cases, be 
metastable [14]. 

In the second case one has L > A, so that Id < lh, and 
in films thinner than lh the paraphase is already unstable 
with respect to domains. Certainly, domains would form 
in thicker films too, so there is no chance that the film 
can become homogeneously polarized. Summarizing, we 
see, that the film can transform with increasing thickness 
(i) from paraphase into homogeneous ferrophase at lh < 
I < Id (which in some cases may be metastable), and 
then into domain state at I > Id, or (ii) the paraphase 



goes over directly into domain state at I > Id- 

Now, we would like to see the effect of a finite band gap 
Eg and free carriers in the ferroelectric itself on stability 
of homogeneous state. One sees from Eq. (3) that at 



P > E g /qL , 



(16) 



where E g is the band gap of the FE, the depolarizing 
field in the FE exceeds \E d \ ~ E g /ql, the band bend- 
ing in the film becomes larger than E g /2, and pock- 
ets with degenerate screening carriers form in atomically 
thin layers at the interface [3]. We can rewrite this into 
a more useful form P > e s P a tKd a t(Eg/8irE at ), where 
Pat = q/dat ~ 200^C/cm 2 is the "atomic" polarization, 
E a t = q 2 /dat ~ 10eV the "atomic" energy. For metallic 
contacts (e s <~ 1 — 3, nd at ~ 1), and a typical band gap 
in FE E g = 2 - 3cV the condition is P > 10 - 20^C/cm 2 , 
so it seems as if it can be met in many ferroelectrics of 
interest not very close to the phase transition. However, 
the system will split into domains close to T c and is un- 
likely to reach such a homogeneous state, for the same 
reason as in FE film without electrodes or in a film with 
lightly doped semiconductor electrodes below. 

In the FE slab without electrodes, the screening by 
the carriers is possible only when the depolarizing field 
again exceeds \Ed\ ~ E g /ql, and pockets of screening 
carriers can form in atomically thin layers at the sur- 
face of FE [2]. Before pockets form, the depolarizing 
field is huge, Ed = —inP and it totally suppresses the 
homogeneous polarization. Therefore, the film cannot 
get into monodomain state until the absolute value of 
thermodynamic coercive field becomes > E g /ql. This 
is certainly not possible until such a low temperature 



where \A\ > \Ai 
II ;>- i> ., 



/g3/2£jl/2£ \ 2 /3 

( — " ) • Using the estimate 



P'ati p at ~ q/dl t , we obtain 



3 3/2 E g dat 

2k B T at I 



2/3 



dat 
I 



2/3 



(17) 



where k B is the Boltzmann constant. If this state could 
be reached, the system would experience a discontinu- 
ous transition into a state with large polarization. But 
it could not, since domains appear much closer to T c . 
Indeed, according to Eqs. (11), (17) 



Ad 
At 



1 



fl \ 1/3 

■f) <]. 



(18) 



The result (17) suggests that in typical ferroelectrics thin- 
ner than 1-10 /jm the monodomain state is simply im- 
possible: the equation of state has only a trivial solu- 
tion P = 0. Similar reasoning applies to the case of 
intrinsic semiconductor electrodes [5,4], where E g now 
means the band gap of the semiconductor. These au- 
thors have found that the transition into a homogeneous 
state should become first order but did not realize that it 
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cannot compete with domains. The above results are in 
striking disagreement with claims [8] that carriers can 
stabilize the homogeneous polarization in the few-nm 
thick FE films with semiconductor electrodes. There is 
no comment on a disagreement with earlier calculations 
by Ivanchik and Guro et al. [2] , who obtained the results 
by solving the equation of state, accounting for presence 
of the pockets with degenerate carriers that were disre- 
garded in [8]. 

We have shown that carriers in either ferroelectric it- 
self or semiconductor electrodes are usually insufficient 
to screen the depolarizing field not far from the phase 
transition. Far from the transition other mechanisms of 
screening may operate, e.g. generation of bulk/surface 
charged defects like oxygen vacancies [20] , surface recon- 
struction [21], or local charge disproportionation between 
ion layers [22,20]. Yet, the formation of domains seems 
to be the dominant screening mechanism in ferroelec- 
tric films, at least near the phase transition. Our re- 
sults suggest that nearly unstrained cubic ferroelectrics 
go over into a multidomain state upon cooling to below 
the critical temperature especially easily. In certain cases 
(perhaps, not observed yet) they could transform into a 
monodomain state. In a multidomain state, the electric 
response of the film would be determined by the proper- 
ties of the domain structure, including pinning. In this 
regard, it may be relevant to the present discussion that 
in 12A PbTi0 3 thin films on SrTi0 3 substrate the Ar- 
gonne group has detected domains [23] . In their case the 
top electrode was absent, but we know that the condi- 
tions for the domain formation in non-electrodcd films 
are similar to those in electroded samples (capacitors) 
if one takes into account possible dead layers or finite 
screening length in the electrodes, see above discussion 
of Eq. (11). Theoretically, at least, the ferroelectric do- 
mains may survive down to thicknesses comparable to 
one unit cell. 
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